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CONVERSION  FACTORS,  NON-SI  TO  SI  (METRIC) 
UNITS  OF  MEASUREMENT 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI 
(metric)  units  as  follows: 


Multipl 


cubic  feet 
degrees  (angle) 


inches 

pounds  (mass) 
per  cubic  foot 


0.02831685 

0.01745329 

0.3048 


16.01846 


To  Obtain 


cubic  metres 

radians 

metres 

centimetres 

kilograms  per 
cubic  metre 


NOMOGRAPHIC  RIPRAP  DESIGN 


PART  I:  DEVELOPMENT 


Introduction 

1.  This  report  presents  a  nomographic  solution  procedure  for  the  design 
of  riprap.  The  procedure  incorporates  five  different  riprap  design  methods 
and  adjusts  riprap  size  for  (a)  specified  safety  factor,  (b)  bend  placement, 
(c)  side  slope  placement,  (d)  a  range  of  riprap  specific  weights,  (e)  stone 
shape,  (f)  use  of  hydraulic  radius  or  local  depth,  and  (g)  specified  roughness 
size.  Examples  are  presented  to  illustrate  the  use  of  the  nomographs. 

Background 

2.  The  development  of  the  five  methods  of  riprap  design  is  presented  in 
detail  elsewhere  (Reese  1984)  and  will  only  be  summarized  here. 

Isbash  method 

3.  By  summing  the  forces  on  a  stone,  Isbash  (1932,  1936)  arrived  at  an 
equation  which,  presented  In  dimensionless  form,  is 


pu: 


(ys  -  Y)d 


2K, 


50 


(I) 


where 

P 


mass  density  of  water* 
velocity  on  the  stone 
specific  weight  of  stone 
specific  weight  of  water 
median  stone  diameter 
empirical  constant 


4.  In  areas  where  no  velocity  profile  has  been  developed  (e.g.,  below  a 
structure) ,  mean  channel  velocity  V  is  commonly  substituted  for  .  The 
value  of  Kj  is  generally  taken  as  0.86  for  high  turbulence  locations  and  1.2 
for  low  turbulence  locations  (Office,  Chief  of  Engineers  (OCE) ,  US  Army, 

1970) . 

Tractive  force-logarithmic  profile 

5.  The  average  boundary  shear  stress  can  be  approximated  by  the  follow¬ 
ing  equation  (OCE  1970) 


„„  ,  ,  12 . 2pD' 

32 ’ 6  lo8l0  “ k 
V  s 


where 

2 

tq  *  average  boundary  shear,  lb /ft 
p  *  multiple  of  depth  such  that  hydraulic  radius  R  =  pD' 

D'  *  depth  to  a  fixed  datum,  ft 

k  ■  equivalent  sand  grain  roughness  of  boundary,  ft 
s 

6.  This  profile  was  originally  derived  using  hydraulic  radius  R 
instead  of  depth  D'  (Keulegan  1938).  Hydraulic  radius,  when  used  with  mean 
channel  velocity,  accounts  somewhat  for  the  channel  shape  effects.  Hydraulic 
radius  is  replaced  by  a  multiple  of  depth  (i.e.,  R  *  pD',  see  Plate  11)  in 
this  and  later  equations.  The  value  of  p  is  taken  equal  to  1  when  local 
mean  depth  is  used  with  the  corresponding  mean  velocity  in  the  vertical  v  to 
approximate  shear  at  a  point  (t  )  .  For  all  irregularly  shaped  cross  sec¬ 
tions  it  is  recommended  that  local  mean  velocity  above  the  stone  be  estimated 
and  used  instead  of  mean  channel  velocity.  See  Table  1  for  further  details. 

7.  Depth  D’  is  measured  from  the  free  surface  down  to  the  assumed 

origin  of  the  logarithmic  profile.  The  original  tractive  force-logarithmic 

profile  (OTFL)  assumed  depth  to  top  of  stone  D  for  D'  and  d^^  for  kg 

(OCE  1970).  Evidence  exists  that  the  actual  origin  of  the  profile  may  be 

about  d  below  the  top  of  the  rock  layer  and  that  k  is  somewhat  greater 

than  dj..  .  Using  known  riprap  roughness  measurements,  k  can  be  shown  to 
jU  s 

be  about  3d^^  (although  the  value  probably  varies  with  depth,  gradation,  and 
shape)  (Reese  1984,  Hey  1979).  U  will  be  used  for  D'  throughout  the  re¬ 
mainder  of  this  report. 


Modified  tractive 
force-logarithmic  profile 


8.  A  modified  profile  (MTFL)  can  be  developed  for  calculation  purposes 
as 
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12. 2p 
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where  m  is  an  empirical  multiplier  as 


k 
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md50 


Again,  for  local  mean  velocity  p  ■  1  and  local  depth  is  used. 

9.  The  critical  tractive  force  is 


(4) 


Tc  *  VYs  '  Y)d50  C5) 

where  K^.  is  Shield's  constant.  The  value  of  may  range  from  about  0.027 

to  0.06  depending  on  flow  conditions  and  definition  of  incipient  motion. 

10.  At  incipient  motion  T0  3  Tc  •  The  OTFL  (Equation  6)  and  MTFL 
(Equation  7)  methods  can  then  be  derived  as 


and 
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Tractive  force-power  profile 

11.  A  power  formula  for  tq  can  be  formed  which  closely  approximates 
the  logarithmic  equation  over  the  range  3  £  R/kg  -  500.  It  was  shown  (Reese 
1984)  that  this  equation  actually  fit  available  data  better  than  the 
logarithmic  profile  and  could  be  solved  directly,  whereas  the  logarithmic 


I 


equation  requires  an  iterative  procedure.  The  power  form  of  is 


63.39 


12.  Combined  with  Equation  5,  the  final  dimensionless  form  is 


(Y_  -  Y) d._  =  63-39KT  mde 


Froude  number  method 


13.  The  Froude  number  method  is  commonly  expressed  as 


f  •  V3 


where 

=  empirical  constant 

F  =  flow  Froude  number  (V/VglT  where  g  =  acceleration  of  gravity) 

14.  When  this  method  is  used,  great  reliance  is  placed  on  model  studies 
and  field  experience  for  the  determination  of  Kp  .  It  is  often  used  in  areas 
of  high  or  uncertain  turbulence  and  in  areas  where  flow  configuration  is  con¬ 
trolled  by  a  rapidly  changing  boundary  (e.g.,  below  structures).  Kp  ranges 
from  0.22  for  incipient  motion  in  flume  studies,  to  0.35  for  channel  design, 
to  as  high  as  7.5  for  low-head  navigation  structures  where  energy  dissipation 
is  poor  (Maynord  1978;  Grace,  Calhoun,  and  Brown  1973).  A  value  of  1.0  is 
recommended  for  structures  with  effective  energy  dissipation. 

15.  Equation  9  can  be  rearranged  as 


( Y  -  Y)63.39KTpI 


This  is  the  same  form  as  Equation  10.  Kp  is  equivalent  to  the  left-hand 
term.  Thus  the  Froude  method  can  be  thought  of  as  a  simplified  tractive 
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K  intrinsically  incorporates  the  p  and  m  adjustments. 

F 

Method  summary 

17.  The  five  equations  for  riprap  design  are  then  (using  V 


Isbash: 


OTFL: 


MTFL : 


_ py: 
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A  table  of  factors  for  converting  non-SI  units  of  measurement  to  SI 
(metric)  units  is  found  on  page  3. 
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PART  II:  SOME  RIPRAP  DESIGN  CONSIDERATIONS 


18.  In  the  design  of  riprap  stone  size,  any  of  three  flow  situations  or 
types  may  predominate.  Depending  on  the  situation,  a  slightly  different  de¬ 
sign  method  is  used.  The  three  types  are  (a)  upstream  roughness  dependence 
(Type  I) ,  (b)  riprap  roughness  dependence  (Type  II) ,  and  (c)  boundary  geometry 
dependence  (Type  III) . 


Type  1 


19.  In  a  Type  I  situation,  flow  moves  from  an  upstream  channel  reach  of 
lesser  roughness  into  the  riprap-protected  reach.  An  example  might  be  transi¬ 
tion  from  a  concrete-lined  section  to  a  riprap-protected  section.  If  the  up¬ 
stream  roughness  is  greater  than  that  of  the  riprap.  Type  II  methods  should  be 
used.  It  can  be  shown  that  if  the  riprap  protrudes  above  the  upstream  bed  by 
more  than  a  small  amount  or  if  the  roughness  difference  is  great,  there  will 
be  a  transition  section  in  which  the  forces  attempting  to  dislodge  the  stone 
will  be  greater  than  further  downstream  where  the  velocity  at  the  top  of  the 
level  of  the  rock  has  been  slowed.  In  Type  I  flow,  the  velocity  of  the  up¬ 
stream  reach  should  be  used  to  size  the  stone  in  this  transition  section.  In 
many  cases  the  velocity  entering  the  section  is  sufficiently  close  to  the 
eventual  velocity  within  the  section  that  Types  I  and  II  riprap  design  method 


will  yield  approximately  equal  d^Q  sizes.  If  the  length  of  the  transition 


section  is  assumed  equal  to  the  distance  required  to  submerse  the  stone  in  a 
turbulent  boundary  layer  initiated  at  the  upstream  edge  of  the  riprap  section, 
the  transition  section  length  can  be  estimated  to  be  75  to  100  times  d-n  • 


20.  Type  II  flow  is  used  if  upstream  and  protected  section  roughness 
are  somewhat  similar  and  riprap  does  not  greatly  protrude  Into  the  flow  (thus 
creating  a  smooth  transition) ,  or  for  protected  sections  downstream  from  a 
transitional  section  when  the  velocity  in  the  channel  and  the  size  of  stone 
are  interdependent.  In  this  case  an  iterative  solution  is  required.  A  size 
Is  first  chosen  and  velocity  calculated.  Then  this  velocity  is  used  to  deter¬ 
mine  a  stable  riprap  size.  The  initial  and  calculated  sizes  are  compared  and 


the  procedure  repeated  until  approximate  agreement  is  reached. 


21.  In  Type  III  flow,  the  velocity  distribution  is  largely  determined 
by  the  boundary  geometry.  Turbulence  caused  by,  for  example,  expansions,  con¬ 
tractions,  drops,  or  sills  predominates.  In  these  cases  the  refinements  of 
velocity  profile  are  overshadowed  by  boundary-generated  macroturbulence.  A 
judicious  choice  of  empirical  constant  is  then  the  most  important  factor  in 
stone  sizing.  The  Froude  or  Isbash  methods  are  usually  used  in  Type  III 
situations  because  of  their  simplicity.  Often  model  studies  are  recommended 
if  the  flow  conditions  are  unique  and  the  importance  of  the  structure  warrants 
the  added  cost. 

Strip  Sections 

22.  Often  the  channel  perimeter  is  not  uniformly  protected.  The  banks 
may  be  protected  and  the  bed  left  in  its  natural  or  excavated  state.  If  a 
large  side  slope  is  to  be  protected,  the  riprap  size  may  be  decreased  in 
longitudinal  strips  up  the  bank  as  a  cost-saving  measure. 

23.  Although  the  actual  scour  mechanism  is  not  well  understood  In  these 
cases,  the  toe  of  the  bank  or  the  lower  portion  of  the  riprap  strip  is  usually 
most  vulnerable  to  scour  or  stone  displacement.  Some  methods  for  toe  protec¬ 
tion  are  discussed  in  Engineer  Manual  1110-2-1601  (OCE  1970).  Actual  required 
toe  depth  Is  quite  variable  and  depends  on  stream  type,  flow  angle,  bed  and 
bank  material,  and  history.  It  is  sometimes  calculated  to  be  the  sum  of  gen¬ 
eral  and  local  scour  plus  one-half  the  anticipated  bed  form  height. 

24.  The  normal  method  of  approach  (short  of  two-dimensional  mathemati¬ 
cal  modeling  but  beyond  simply  using  mean  channel  velocity)  when  roughness 
varies  in  longitudinal  strips  is  first  to  calculate  the  mean  velocity  in  each 
strip  (e.g.,  using  a  suitable  backwater  program  for  gradually  varied  flow  or 
the  alpha  method  (OCE  1970)  for  uniform  flow) .  It  is  then  assumed  that  the 
faster  velocity  of  an  adjacent  strip  impinges  on  a  section  of  slower  velocity 
at  unknown  and  varying  points.  For  example,  the  faster  velocity  in  the  center 
portion  of  a  channel  with  protected  banks  is  used  to  size  the  riprap  on  the 
bank,  or  the  faster  velocity  on  a  lower  portion  of  a  bank  is  used  to  size  the 
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riprap  on  an  adjacent  upper  portion  of  the  bank.  These  methods  are  illus¬ 
trated  in  Table  1. 

25.  The  built-in  safety  factor  incurred  by  using  this  method  will,  for 
steeper  slopes  and  shallower  flows,  give  unrealistically  large  riprap  sizes. 
For  certain  combinations  of  variables  no  solution  can  be  found  for  the 
logarithmic  methods.  It  is  probably  best,  in  these  situations,  to  reconsider 
whether  riprap  is  an  appropriate  form  of  bed  or  bank  protection  and,  if  it  is, 
to  choose  a  size  based  on  experience  or  velocity  measurement. 


PART  Ills  THE  NOMOGRAPHIC  METHOD 


Basis  of  the  Nomographic  Method 


26.  The  Nomographic  Method  is  an  attempt  to  provide  a  simple-to-use  yet 
accurate  method  for  riprap  sizing.  It  has  an  additional  advantage  in  that  it 
also  provides  a  design  engineer  with  a  visual  "feel"  for  the  relative  sensi¬ 
tivity  of  various  assumptions  and  placement  situations  used  in  design. 

27.  The  basic  concept  of  this  method  is  illustrated  by  Plate  1.  For  a 
given  set  of  conditions,  a  velocity  versus  depth  plot  can  be  drawn  using  a 
family  of  d^^  curves.  The  designer  simply  finds  the  intersection  of  known 
velocity  and  depth  and  reads  the  appropriate  d,.^  size  directly.  Plate  1 
uses  the  OTFL  equation  (Equation  6)  for  a  stated  set  of  "base  conditions." 

The  problem  is  that  there  is  an  infinite  number  of  combinations  of  conditions 
each  of  which  would  generate  a  different  family  of  curves.  For  example,  if 
the  riprap  is  placed  on  a  side  slope,  the  family  of  curves  would  shift  to  the 
left,  requiring  a  larger  d,.g  size  for  the  same  velocity  and  depth. 

28.  Rather  than  shifting  the  d^Q  curves  to  the  left  (or  right)  when 
departures  from  base  conditions  are  encountered,  in  the  nomographic  method, 
the  velocity  is  shifted.  The  question  can  be  asked,  "What  velocity  would  re¬ 
sult  in  the  same  d^g  size  using  the  base  condition  family  of  curves  (Plate  1 
here)  as  would  be  obtained  by  plotting  a  new  family  of  curves  for  the  new  set 
of  conditions?" 

29.  A  departure  from  any  of  the  base  conditions  (Plates  1-5)  is  ac¬ 
counted  for  by  an  adjustment  to  the  real-world  velocity,  creating  a  fictitious 
velocity.  That  is,  any  factor  which  increases  or  decreases  the  ability  of  a 
stone  to  remain  stationary  can  be  reflected  by  a  corresponding  decrease  or 
increase  in  velocity.  This  fictitious  new  velocity  is  called  the  "effective 
velocity."  Thus,  through  coefficients  obtained  from  Plates  6-13,  the  user 
moves  either  to  the  right  or  left  of  the  actual  velocity  on  the  velocity  scale 
and  then  moves  upward  to  the  calculated  depth  and  reads  the  d^g  size. 

30.  The  base  conditions  (Plates  1-5)  can  be  adjusted  for  riprap  shape, 
side  slope  or  bend  placement,  specific  weight  changes,  use  of  different  con¬ 
stant  values  or  safety  factors,  use  of  hydraulic  radius  (channel  shape),  and 
use  of  different  equivalent  sand  grain  roughnesses.  Examples  in  Table  l 
illustrate  adjustments  to  the  base  conditions.  The  following  base  conditions 


Rk 


ora 


/.  -* 


are  used  in  Plates  1-5:  moderately  angular  riprap  placed  on  the  bed  of  a 
straight  channel,  Yg  -  165  lb/cu  ft  ,  kg  =  d^  ,  D  is  used  instead  of  R  , 
and  base  values  for  the  constants  are  used  (R^  =  0.04  ,  R^  -  0.22  ,  R 
=  1.2  ). 

31.  The  tractive  stress  displacement  methods  depend  on  an  accurate 
description  of  the  vertical  velocity  profile  to  be  able  to  transpose  a  known 
velocity  to  a  bed  shear  stress.  When  the  bed  roughness  elements  are  large  in 
relation  to  the  flow  depth  (between  1/10  to  1/3),  the  methods  begin  to  lose 
accuracy  due  to  roughness  shape  and  free  surface  effects.  These  effects  tend 
to  increase  the  friction  factor  (slow  the  velocity)  and  thus  lead  to  a  more 
conservative  design.  The  dashed  lines  on  Plates  1-5  indicate  a  D/d,.Q  va^-ue 
of  3  (roughness  elements  *  1/3  depth).  These  charts  should  not  be  used  for 
design  below  this  value  as  variation  from  profile  assumptions  reaches  unac¬ 
ceptable  levels. 


Derivation 


32.  At  incipient  motion  T0  =  Tc  •  Several  factors  may  increase  actual 

shear  (t  )  or  decrease  critical  shear  (t  ). 
o  c 

Bend  placement 

33.  The  maximum  shear  in  a  bend  (OCE  1970)  can  be  expressed  as 


•  where 

t.  *  shear  in  a  bend 

P  D 

T  -  top  width  of  channel  entering  bend 

R  ■  centerline  radius  of  curvature  of  bend 
c 

or 


(15) 


(16) 


where  ,  the  coefficient  for  bend  shear,  is  the  term  in  brackets  in 
Equation  15. 


34.  The  required  critical  shear  for  side  slope  placement  (OCE  1970)  is 
expressed  as 


sin^  $ 


.  2  „ 
sin  0 


where 


t  *  x  K 
s  c  s 


rg  ■  side  slope  critical  shear 
=  side  slope  angle 
0  *  riprap  angle  of  repose 

t  *  coefficient  of  side  slope  shear  (bracketed  term  in  Equation  17) 


Safety  factor 


35.  The  safety  factor  is  expressed  as 


ToKSF 


where 


Tgp  *  actual  shear  with  a  safety  factor  applied 


Kgp  -  safety  factor 

A  factor  of  1.5  had  been  recommended  for  moderate  nonuniformity  in  plan, 
profile,  or  roughness  (OCE  1971). 

36.  All  of  these  factors  can  be  applied  independently  (although  the 
bend  correction  equation  incorporates  about  a  1.5  safety  factor  with  it  when 
compared  with  straight  channel  design)  at  the  discretion  of  the  designer. 
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37 .  Remembering  that  k 


and  R  *  pD  and  using  to  repre¬ 


sent  a  specific  weight  of  stone  other  than  165  lb/cu  ft  (base  condition) , 
Equations  5  and  8  are  substituted  into  21  for  and  ,  respectively, 

rearranged  to  give 


[-63.39^;  -  y)d50Kl  /2  /  \ 

L  KBP  J  \d50/ 


where  is  the  use  of  a  constant  other  than  the  base  constant.  K^,  com¬ 
bines  Kt  and  KgF  such  that  K.J,  -  •  Tllus  a  choice  of  K^,  of  0.027 

is  equivalent  to  a  shear  safety  factor  of  1.5  (0.027  *  0.04/1.5)  and  so  on. 

38.  The  variable  V  in  Equation  22  is  the  actual  velocity.  The 
effective  velocity  is  that  which  will  give  the  same  d^^  for  the  same  depth 
but  for  base  conditions.  It  is  expressed  as 


1  1/2  /  ■ 

]  te 


'39V*,  '  *>*50 


where  the  base  value  for  K_  is  0.04  and  y  is  165  lb/cu  ft. 

l  s 

39.  Dividing  Equation  23  by  22  gives 


Vf  •  /o.O«V/Z  r  102,6  1  1/2  (h\  1/2  (m\ 

v  *  vv)  L<*;  -  *> J  w 


V  __  -  C  C  C.  C  C  C  V 
eff  t  g  b  s  m  p 
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Thus  for  any  given  set  of  conditions  the  user  determines  the  correction 
coefficients  and  multiplies  them  times  the  given  velocity  to  obtain  an 


effective  velocity.  He  then  reads  the  d^g  size  directly  from  the 


intersection  point  of  effective  velocity  and  depth. 


Froude  and  Isbash  Methods 


40. 


Because  of  the  nature  of  the  Froude  and  Isbash  methods,  the  C 


and 


corrections  are  not  applicable  to  them.  Side  slope  and  bend 


corrections  are  not  applicable  to  the  Isbash  method  in  its  original 
derivation. 


Logarithmic  Profile 


41.  The  logarithmic  profile  can  be  arranged  similar  to  the  power 
profile  giving 


17 


h 


Veff  '  5’75 


V^s  -  Y)df 


!J 


Dividing  Equation  27  by  Equation  26  gives 


f(102.6  )11/2  (h\U2 

J  \h) 
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V  ,  «  C  C  C.C  SV 

eff  t  g  b  s 

where  ,  and  Cg  are  as  defined  following  Equation  25  and 
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C  carries  with  it  the  correction  for  both  hydraulic  radius  and  equivalent 
roughness  and  was  arranged  in  this  way  to  provide  multiplicative  rather  than 
additive  corrections.  For  local  velocity  use  set  R  =*  D  (the  extreme  top  of 
the  chart) . 


Summary 


42.  The  various  coefficients  used  to  adjust  the  actual  velocity  can  now 


*Vr* 


be  summarized.  These  equations  may  be  used  in  lieu  of  the  plates  if  greater 
accuracy  is  desired. 

a.  Stone  shape  and  side  slope  placement  C  corrections  are  found 
in  Plates  6  and  7,  respectively,  and  expressed  as 


2 

sin  © 


Plate  6  was  developed  from  work  by  Simons  (1957)  and  is  used  to 
estimate  angle  of  repose. 

b.  The  bend  correction,  found  in  Plate  8,  is  expressed  as 


S  •  »•”  r 


c.  The  specific  weight  correction  is  found  from  Plate  9  and 
expressed  as 


102.6 

(Y*  *  62.4)  _ 


d.  The  constant  correction  coefficient  for  each  of  the  four 

methods  is  found  in  Plate  10.  For  the  power  and  logarithmic 
profiles  it  is  given  as 


For  the  Isbash  method  it  is 


C  *  14 
i 


where  K|  is  the  chosen  Isbash  constant.  For  the  Froude 
method  it  is 


where 


'  is  the  chosen  Froude  constant, 
r 

e.  The  channel  shape  correction  factor  for  the  power  profile  is 
found  in  Plate  11  for  a  known  top  width,  depth,  and  side  slope. 
It  is  given  as 


C 

P 


(37) 


£.  The  correction  for  both  shape  of  channel  and  equivalent  rough¬ 
ness  for  the  logarithmic  profiles  is  given  in  Plate  12  and  as 


C 


log10 

log10 


(38) 


The  depth  used  here  depends  on  which  logarithmic  method  will  be 
used.  The  OTFL  method  uses  depth  to  top  of  rock  while  the  MTFL 
method  uses  depth  to  d^^  below  top  of  rock.  Therefore  use 

D/d^0  +1  in  place  of  D/d^Q  in  both  numerator  and  denomi¬ 
nator  of  Equation  38  for  the  MTFL  method. 

£.  Plate  13  gives  the  power  profile  correction  for  equivalent 
roughness  as 


C 

m 


m 


1/6 


(39) 


Safety  Factors 


43.  Safety  factors  (OCE  1971)  for  the  logarithmic  and  power  profiles 
are  found  from 


SF 


T 


0.04 


(40) 


where  SF^  is  the  shear  safety  factor.  For  the  Froude  method,  assuming  0.22 
to  be  the  ideal  situation  incipient  motion  constant,  it  is  found  that 


2/3 


These  are  shear  safety  factors.  For  example:  to  provide  the  recommended  1.5 
safety  factor  on  the  calculated  actual  shear,  must  be  0.027  and  K^,  must 


be  0.404. 


Shear  Analysis 


44.  Insight  can  be  gained  by  returning  to  the  original  premise  for 

shear-induced  riprap  instability  of  paragraph  10.  At  incipient  motion  the 

shear  acting  on  a  stone  equals  the  critical  shear  for  movement  or  x  -  x  . 

o  c 

Equation  20  demonstrated  in  a  preliminary  way  the  relationship  between  the 
actual  and  critical  shears  xq  and  x£  ,  and  the  various  factors  affecting 
them.  By  combining  Equations  31-38  (the  factors  read  directly  from  the 
nomograph)  with  the  definitions  for  the  shears,  a  nomographic  shear  analysis 
can  be  performed,  identical  to  the  method  given  in  EM  1110-2-1601  (OCE  1970). 

45.  For  the  OTFL  method  the  proper  equation  is 


<VtC>2 '  4’10d50  fe) 2 

10  dcrt 


32.6  log 


Here,  and  in  the  two  equations  following,  the  form  of  the  equation  is 


The  left  side  is  the  actual  shear  adjusted  for  channel  shape,  bend,  equivalent 
roughness,  and  safety  factor.  The  right  side  of  the  equation  is  the  critical 
shear  adjusted  for  specific  weight  and  side  slope  placement. 

46.  The  equation  for  the  MTFL  method  is 


32.6  log1Q  12.2 
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47.  The  equation  for  the  power  profile  method  is 
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Computer  Programs 

48.  A  computer  program  which  allows  for  on-line,  interactive  riprap 
design  using  the  methods  presented  in  this  paper  is  available  at  all  US  Army 
Corps  of  Engineers  design  offices  through  the  Conversationally-Oriented  Real- 
Time  Programming  System  (CORPS).  It  is  entitled  H701 1-Riprap  Design  by  Four 
Methods.  For  complete  documentation  and  information,  contact  the  Engineer 
Computer  Programs  Library  (WESIM-SC)  at  FTS  542-2581. 
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Table  1 
Examples 


Example  1:  Type  I  Flow — Bend  Placement 


Given: 


R  =  200  ft 
c 

T  *  50  ft 


D  *  8  ft 


V  *  11  fps 


y  =  175  lb/cu  ft 
s 


Side  Slope  *  2H:1V 

Moderately  rounded  stone,  trapezoidal  channel 
Use  TFP  method  and  mean  channel  velocity  (upstream) 

Step  1.  Side  Slope  Correction 

Assume  d^Q  *  15  in.  $  *  26.6  deg  *  0.027  (1.5  safety  factor) 

From  Plate  6:  0  =  40.5  deg 

From  Plate  7:  C  =1.17 
s 

Step  2.  Bend  Correction 

T/R  =  50/200  =  0.25 
c 

From  Plate  8:  C,  =  1.25 

D 

Step  3.  Specific  Weight  Correction 

From  Plate  9:  C  *  0.95 

g 

Step  4.  Power  Profile  Channel  Shape  Correction 

T/D  =  50/8  «  6.25 

From  Plate  11:  C  =  1.07 
P 

Step  5.  Equivalent  Sand  Grain  Roughness  Correction,  Power  Profile 
m  *  2 

From  Plate  13:  C  =1.12 
m 

Step  6.  Safety  Factor 
K^,  =  0.027 


(Continued) 
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Table  1  (Continued) 


From  Plate  10:  Cfc  -  1.22 


Step  7.  Power  Profile 

V  ,,  =  C  C.C  C  C  C  V 
eff  s  b  g  p  m  t 

=  (1.17) (1.25) (0.95) (1.07) (1.12)  (1.22)  (11) 
=  22.34  fps 

From  Plate  4:  d,.g  *  31  in. 

Step  8.  Check  Side  Slope  Correction 


d5Q  =  31  in. 


$  *  26.6  deg 


k;  -  0. 


From  Plate  6:  0  *  41.0  deg 


From  Plate  7:  C  =1.17 
s 


The  side  slope  correction  does  not  change  so  the  rock  size  dqn  =  31  in.  is 
correct. 

Example  2:  Type  I  Flow — Transition  Design 


Given:  Q  ■  4,000  cfs 


Bottom  Slope  =  0.00082 


A.  Concrete  Section 


Manning's  n  ■  0.012 


D  =  6  ft 


B.  Transition  Section 


Bottom  width  «  50  ft  “  165  lb/cu  ft 

2H:1V  Side  Slopes 


V  (mean)  ■  10.75  fps 

V  (center)  =  11.69  fps 

V  (sides)  =  6.84  fps 


Assume  a  fairly  smooth  transition 


Power  method 


*  0.027  (1.5  safety  factor) 


Step  1.  Side  Slope  Correction 
Assuming  d^Q  =  24  in. 


(Continued) 
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From  Plate  6:  0  =  41.5  deg  (moderately  angular  stone) 

From  Plate  7:  C  =1.16 
s 


Step  2.  Safety  Factor 


From  Plate  10:  =  1.22 


Step  3.  k  *  3dco 
r  s  50 


From  Plate  13:  C  *  1.20 
m 

Step  4.  V  =  V(C  )(C  )(C  )  =  (11.69)  (1. 16) (1 .22) (1. 20)  =  19.85  fps 
et  t  s  t  m 

From  Plate  4:  d^^  *  25  in. 

Step  5.  Assume  approximate  transition  zone  distance  of  75d,.Q  =  150  ft  . 
C.  Downstream  Section  Size  Check 


Assume  n  *  0.030  (conservative) 
Then  V  (center)  =  6.54  fps 


V  (mean)  =  5.77  fps 


D  =  9.93  ft 


3.5  in. 


D.  The  designer  now  has  the  choice  of  how  to  step  down  the  size  from 
24  in.  to  an  unprotected  section.  The  24-in.  size  will  be  carried  100  ft 
downstream  with  a  12-in.  stepdown  for  an  additional  50  ft  with  tie-ins. 

Example  3:  Type  II  Flow — Side  Slope  Placement 


Given:  Q  =  8,000  cfs 

2H:1V  side  slopes 


Bottom  width  =  70  ft 


Bottom  slope  *  0.0040 


Bottom  roughness:  Manning's  n  =  0.028 

=  165  lb/cu  ft  (moderately  angular  stone) 


Assume  uniform  flow 


1.5  safety  factor 


k  =  2d,0 
s  50 


(Continued) 
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Table  1  (Continued) 


Size  stone  for  bank  protection  using  power  profile  method  and 
sectional  velocities. 

Step  1.  Determine  applicable  velocities  and  depth.  Using  CORPS  program  H7012 
-  The  Alpha  Method: 


Assuming  Manning's  n  =  0.04dt 


(d^Q  in  feet)  and  d5Q  =  12  in. 


Manning's  n  *  0.04  on  protected  slopes 

Center  velocity  *  13.28  fps 

Depth  =  7.88  ft 

Side  slope  velocity  =  5.44  fps 

Step  2.  Side  slope  correction 

From  Plate  6:  Repose  angle  ■  41  deg 

From  Plate  7:  C  -  1.17 
s 

Step  3.  Equivalent  roughness  correction 

From  Plate  13  (m  =  2) :  C  =1.12 

m 

Step  4.  Safety  factor  correction 

From  Equation  40  for  SF^  =  1.5:  K^,  =  0.027 

From  Plate  10:  Ct  ■  1.22 

Step  5.  V  =  (1.17)  (1.12) (1.22) (13.28)  =  21.32  fps 
From  Plate  4:  d^g  ■  26  in. 

Step  6.  Recalculate  downstream  velocity  using  larger  d^g  size 
Manning's  n  =  0.04(26/12)^^  =  0.045 
Using  H7012 

Velocity  center  =  13.33  fps 
Velocity  bank  =  4.85  fps 
Depth  =  7.92  ft 

(Continued) 
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Step  7.  Recalculate  size 

V  *  (1.17)(1.12)(1. 22) (13.33)  -  21.31  fps 
et  t 

From  Plate  4:  *  27  in.  (OK) 

Example  4:  Type  III  Flow — Basin  Design 
Flow  below  drop  structure.  Use  Isbash  and  Froude  methods. 

V  (over  end  sill)  =  12.4  fps 
D  =  11  ft  y  -  185  lb/cu  ft 


3H:1V  side  slopes  with  moderately  angular  stone 
and  K^,  are  selected  based  on  experience 
Step  1.  Specific  weight  correction 


From  Plate  9:  C  *  0.91 

g 


Step  2.  Side  slope  correction  (Froude  only) 
From  Plate  6:  0  *■  41.5  deg 


From  Plate  7:  C  *  1.07 
s 


Step  3.  From  Plate  10:  =*  1.39  for  »  0.86  (for  high-turbulence 

location) 

C,  *  1.65  for  Kj,  *  1.00  (for  structures  with 


effective  energy  dissipation) 
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BASE  VALUES 
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DASHED  LINE  DUE  TO  FREE  SURFACE 
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CRITICAL  MEAN  VELOCITY 
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PLATE  4 


1.6 


SPECIFIC  WEIGHT 
CORRECTION  COEFFICIENT 


CORRECTION  COEFFICIENT 
LOGARITHMIC  PROFILE 


7/y/ 


OBTAIN  EFFECTIVE  VELOCITY  CORRECTION  COEFFICIENT 

POWER  PROFILE 


APPENDIX  A:  NOTATION 


Logarithmic  profile  modification  coefficient  for  hydraulic  radius 
and  equivalent  roughness 

Coefficient  for  velocity  modification  due  to  a  Shields  constant 
using  a  bend  shear 

Correction  coefficient  for  Froude  method 

Coefficient  for  velocity  modification  due  to  a  Shields  constant 
using  a  specific  weight 

Correction  coefficient  for  Isbash  method 

Coefficient  for  velocity  modification  due  to  a  Shields  constant 
using  a  equivalent  roughness 

Coefficient  for  velocity  modification  due  to  a  Shields  constant 
using  a  hydraulic  radius 

Coefficient  for  velocity  modification  due  to  a  Shields  constant 
using  a  side  slope  shear 

Coefficient  for  velocity  modification  due  to  a  Shields  constant 
using  a  safety  factor 

Median  stone  diameter 

Depth  to  top  of  stone 

Depth  to  a  fixed  datum,  ft 

Flow  Froude  number  (V/VgF) 

Acceleration  of  gravity 

Equivalent  sand  grain  roughness  of  boundary,  ft 
Coefficient  of  side  slope  shear 
Coefficient  for  bend  shear 
Empirical  constant 


Chosen  Froude  constant 


Empirical  constant 


Chosen  Isbash  constant 


,  /.  w*4  ■c,  A  ■*,  •*.  »"« •  *. ,*11* »  V .  ^ hj' ^ ^ k  ^ 
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Safety  factor 


i 

*T 

m 

MTFL 

OTFL 


P 

Q 

R 


R 

c 

SFp 

sft 

T 

TFP 

Ub 

v 

V 


eff 

Y 


Cef  f 


Shield's  constant 

Constant  similar  to  except  it  includes  a  safety  factor 

Empirical  multiplier 

Modified  tractive  force-logarithmic  profile 
Original  tractive  force-logarithmic  profile 
Multiple  of  depth  such  that  hydraulic  radius  R  »  pD' 

Discharge 
Hydraulic  radius 

Centerline  radius  of  curvature  of  bend 
Shear  safety  factor  for  Froude  method 

Shear  safety  factor  for  logarithmic  and  power  profile  method 

Top  width  of  channel  entering  bend 

Tractive  force-power  profile 

Velocity  on  the  stone 

Local  mean  velocity  in  the  vertical 

Mean  channel  velocity 

Effective  velocity 

Specific  weight  of  water 

Specific  weight  of  stone 

Specific  weight  of  stone  other  than  base  condition 
Riprap  angle  of  repose 
Mass  density  of  water 

2 

Shear  in  a  bend,  Ib/ft 
Critical  tractive  force 

Critical  shear  adjusted  for  specific  weight  and  side  slope  placement 


Average  boundary  shear,  lb/ft 

Actual  shear  adjusted  for  channel  shape,  bend,  equivalent  roughness 
and  safety  factor 

2 

Side  slope  critical  shear,  lb/ft 
Actual  shear  with  a  safety  factor  applied 
Side  slope  angle 
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